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Abstract 

The exp(-0(/7))-expansion method is an ascending method for obtaining exact and solitary wave solutions for 
nonlinear evolution equations. In this article, we implement the exp(-0(/7))-expansion method to build solitary 
wave solutions to the fourth order Boussinesq equation. The procedure is simple, direct and useful with the help of 
computer algebra. By using this method, we obtain solitary wave solutions in terms of the hyperbolic functions, the 
trigonometric functions and elementary functions. The results show that the exp(-0(/7))-expansion method is 
straightforward and effective mathematical tool for the treatment of nonlinear evolution equations in mathematical 
physics and engineering. 
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Background 

The world around us is inherently nonlinear (He 2009) 
and nonlinear evolution equations (NLEEs) are widely 
used as models to describe complex physical phenomena 
in various fields of science and engineering, especially in 
solid-state physics, plasma physics, fluid mechanics, 
biology etc. One of the fundamental problems for these 
models is to obtain their travelling wave solutions as 
well as solitary wave solutions. In particular, various 
methods have been utilized to explore different kinds of 
solutions of physical problems described by nonlinear 
evolution equations. In the numerical methods, stability 
and convergence should be considered, so as to avoid 
divergence or inappropriate results. However, in recent 
times, a variety of analytical and semi-analytical methods 
have been developed and use for solving NLEEs, for 
instance, the inverse scattering transform (Ablowitz and 
Clarkson 1991), the complex hyperbolic function method 
(Chow 1995; Zayed et al. 2006), the rank analysis 
method (Feng 2000), the ansatz method (Hu 2001a, b), 
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the (GVG)-expansion method (Wang et al. 2008; Bekir 
2008; Neyrame et al. 2012; Akbar et al. 2012; Alam and 
Akbar 2013; Alam et al. 2014), the Exp-functions 
method (He and Wu 2006), the modified simple equa- 
tion method (Jawad et al. 2010; Khan et al. 2013), the 
Jacobi elliptic function method (Chen and Wang 2005; 
Liu 2005), the Adomian decomposition method (Adomian 
1994; Wazwaz 2002), the homogeneous balance method 
(Wang 1995; Zayed et al. 2004), the F-expansion method 
(Wang and Zhou 2003; Wang and Li 2005), the Backlund 
transformation method (Miura 1978), the Darboux 
transformation method (Matveev and Salle 1991), the 
homotopy perturbation method (Mohyud-Din 2007; 
Mohyud-Din and Noor 2009), the generalized Riccati 
equation method (Yan and Zhang 2001), the tanh-function 
method (Wazwaz 2005), the Hirotas bilinear method 
(Hirota 2004), the auxiliary equation method (Sirendaoreji 
2007), the exp(-0(^))-expansion method (Khan and 
Akbar 2013) etc. 

The objective of this article is to implement the potential 
exp(-0(^))- expansion method to search solitary wave 
solutions for nonlinear evolution equations via the fourth 
order Boussinesq equation. In former literature, the 
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solitary wave solutions to the Boussinesq equation 
have not been studied by this method. 

The article is organized as follows: In Methodology 
Section, we give the description of the exp(-0(^))- ex- 
pansion method. We apply this method to the fourth 
order Boussinesq equation in the Application Section. 2D 
and 3D graphs are given in the Graphical Representations 
of the Solutions Section. Finally, in Conclusion Section, 
we draw our conclusions. 

Methodology 

Let us consider the nonlinear evolution equation in 
the form 



F(li^ U x ^ U xx ^ Utti Ut x ^ ...), 



(1) 



where u = u(x, t) is an unknown function, F is a polyno- 
mial in u(x, t) and its derivatives in which highest order 
derivatives and nonlinear terms are involved and the 
subscripts indicate partial derivatives. In order to in- 
vestigate solitary wave solutions of (1) by using the exp 
(-0(^))-expansion method, we have to perform the 
following important steps: 

Step 1. We combine the real variables x and t by a 
compound variable rj 



u{x,t) = u(rj), tj = x±Vt, 



(2) 



where V is the celerity of the traveling wave. By means 
of traveling wave transformation (2), Eq. (1) switch into 
an ordinary differential equation (ODE) for u - u(rj): 



(3) 



where H is a polynomial of u and its derivatives and 
the superscripts refer to the ordinary derivatives with 
respect to tj. 

Step 2. Assume the traveling wave solution of (3) can 
be articulated as follows: 



"(v) = £^i(exp(-cD(/7))) 



(4) 



A t (0 < i < N) are constants to be determined, such 
that A N * 0 and O = O(^) satisfies the following auxiliary 
equation: 



Q>\rj) = exp(-(D(^)) +^exp(0(^)) + A, 



(5) 



Depending on the parameters involved, Eq. (5) has the 
subsequent solutions: 
When pi * 0, and A 2 - 4^ > 0, 



0>{rj) = In 



^-yV-^tanh {rj + E fj _ X 



\ 



2{t 



(6) 



When p * 0, and A 2 - 4^ < 0 ; 
0>{rj) = In 



/■/(^)tanfv^ 



V 



2fi 



(V + E) )-A^ 



When fi = 0, A * o, and A 2 - 4// > 0, 
A 



<D(//) = - In 



s exp(A(>7 + £))-l 
When ft * 0, A * 0, and A 2 - 4// = 0 
2(A(^ + £)+2) 



<D(//) = In 



A 2 (rj + E) 



When fi = 0, A = 0, and A 2 - 4fi = 0, 
<D(//) = ln(// + £) 



(7) 



(8) 



(9) 



(10) 



Step 3. The positive integer N can be determined 
by considering the balance between the highest order 
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derivatives and the nonlinear terms of the highest order 
appearing in (3). 

Step 4. We substitute Eq. (4) into Eq. (3) and then we 
take into consideration the function exp(-0(^)). In con- 
sequence of this substitution, we obtain a polynomial in 
exp(-0(^)). We collect all the coefficients of identical 
power of exp(-0(^)) and equalize to zero delivers a sys- 
tem of algebraic equations whichever can be solved to 

find A N , ,KA,^. The values of A N) ,KA,^ along 

with general solutions of Eq. (5) complete the determin- 
ation of the solution of Eq. (1). 

Application 

In this section, we will use the exp(-0(^))-expansion 
method to construct the exact solutions and then the 
solitary wave solutions to the fourth order Boussinesq 
equation. Let us consider the equation 



u tt~ u x 



-3(u z 



0. 



(11) 



The above model (11) was introduced by Boussinesq 
to illustrate the propagation of long waves in shallow 
water (Lai et al. 2008), where u(x, t) is the elevation of 
the free surface of the fluid, where the subscripts denot- 
ing partial derivatives. The equation also arises in many 
other physical applications, such as, nonlinear lattice 
waves, iron sound waves in plasma, and vibrations in a 
nonlinear string. It was also applied to the study of the 
percolation of water in porous subsurface strata. 

Equation (11) possesses solitary waves, extract from 
traveling wave solutions and Boussinesq was the first who 



gave a scientific explanation of their existence. We utilize 
the traveling wave variable u{rj) = u(x, t), rj = x - Vt and this 
operation changes (11) to the following ODE: 

V 2 u-u"-u^-3{u 2 )" = 0. (12) 

Integrating Eq. (12) twice with respect to rj yields: 

(V 2 -l) u-u -3u 2 + C = 0, (13) 

where C is an integration constant to be determined. 

Balancing the highest order nonlinear term u 1 and linear 
term of the highest order u appearing in (13), yields N = 2. 
Therefore, the solution of Eq. (13) takes the form 



u(rj) = A 0 + Ai exp(-® (//)) + A 2 exp(-0(//)) 2 



(14) 



where A 0 , A lf A 2 are arbitrary constants such that A 2 * 0. 

We substitute Eq. (14) into Eq. (13) and taking consid- 
eration Eq. (5), it generates a polynomial and then setting 
the coefficients of exp(-0(^)) to zero, yields 

-6A 2 -3A 2 2 = 0, (17) 

-2Ai-10A 2 A-6AiA 2 = 0, (18) 

-A 2 -4A 2 X 2 -3A l X-3A 1 2 -6A 0 A 2 + V 2 A 2 -8A 2 {t = 0, 

(19) 

C + V 2 A Q -A 0 -2A 2 ii 2 -3A Q 2 -3A l 2 -A l Xii = 0, (20) 
Solutions of Eqs. (17)-(20), yield 
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Figure 4 Sharp bell shaped 2D and 3D-plot of the solitary wave. The graphs have been plotted from u 4 {q) when fj = 0.25, X = 1, E = 1, f = 0 

and A Q = 1. 



C = -4^ 2 -3A 0 2 -2AV-A 2 A 0 -8A 0/ w 7 



V = ±A/ A + 1 + S> + 6A 0 , A a = -2A,Ai 



-2. 



where A, and A 0 are arbitrary constants. 

Substituting the values of V,A 0 ,Ai,A 2 into Eq. (14), 
yields 



u 2 (rj) =A 0 



4X{* 



4//-A 2 tan (// + £) -A 



8^ 2 



4//-A 2 tan 



4//-A 2 



-(// + £) -A 



u(rj) = A 0 -2X exp(-0(//))-2 exp(-0(//)) 2 



(21) 



where 77 = A 2 + 1 + 8^ + 6A 0 £. 

By using the solutions of Eq. (5) into Eq. (21), we obtain 
the succeeding traveling wave solutions of the fourth order 
Boussinesq equation: 

Type 1: When ft * 0, A 2 - 4^ > 0, 



ui(rj) = A 0 +- 



A 2 -4^tanh| 



A -4^ 



-(// + £) +A 



A -4/Uanh 



A 2 -4// 



-(// + £) +A 



where rj = x+yX + 8// + 6A 0 + 1 £ and £ is an arbitrary 
constant. 
Type 2: When ft * 0, A 2 - 4^ < 0, 



where rj = x+yX + 8// + 6A 0 + 1 £ and £ is an arbitrary 
constant. 

Type 3: When \i = 0, A * 0, and A 2 - 4// > 0, 

^ 2A 2 2A 2 

U3{r/) ~ °~ exp(Afa + £))-! " j expUf/y + £))-l| 2 ' 
where rj = x+ ^Jx 2 + 8// + 6A 0 + 1 £ and £ is an arbitrary 
constant. 

Type 4: When * 0, A * 0, and A 2 - 4^ = 0, 



u 4 (rj) =A 0 + 



A 3 (// + £) 1 / X\rj + E) 



X(rj + £) + 2 2 \A(// + £) + 2 



where rj = x+yX + 8// + 6A 0 + 1 £ and £ is an arbitrary 
constant. 

Type 5: When ^ = 0, A = 0, and A 2 - 4^ = 0, 

M5O7) = Ar 



where rj = x+yX + 8// + 6A 0 + 1 £ and £ is an arbitrary 
constant. 
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Figure 5 2D and 3D-plot of the solitary wave. The graphs have been plotted from u 5 {q) when fj = 0, X = 0, E = 1,t = 0 and /\ 0 = 1. 
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Solitary wave solutions represent an important type 
of solutions for nonlinear partial differential equations 
(PDEs) as many nonlinear partial differential equations 
have been found to have a variety of solitary wave solu- 
tions. It is familiar that searching of exact solutions of 
nonlinear partial differential equations plays a significant 
role in the study of nonlinear physical phenomena. Exact 
traveling wave solutions are useful for verifying the accur- 
acy and stability of popular numerical schemes such as the 
finite difference and finite element methods. The solitary 
wave solutions obtained in this article are encouraging, 
applicable, and could be helpful in analyzing long wave 
propagation on the surface of a fluid layer under the ac- 
tion of gravity, iron sound waves in plasma, and vibra- 
tions in a nonlinear string. 

Graphical representation of the solutions 

Solitary waves can be obtained from each traveling wave 
solution by setting particular values to its unknown pa- 
rameters. By adjusting these parameters, one can get 
an internal localized mode. In this section, we have 
presented some graphs of solitary waves constructed 
by taking suitable values of the involved unknown pa- 
rameters to visualize the underlying mechanism to the 
original physical phenomena. Using mathematical soft- 
ware Mathematica, two and three-dimensional plots of 
the obtained solutions have been shown in Figures 1, 2, 
3, 4 and 5. 

The obtained solutions of the fourth order Boussinesq 
equation incorporate four types of explicit solutions 
namely hyperbolic, trigonometric, exponential, and ra- 
tional function solutions. From these explicit results we 
observe that solutions u x (y\) and u 2 {rj) are soliton and 
the rest of the three solutions are cuspon. The above 
solitary wave solutions might be useful in analyzing the 
propagation of long waves in shallow water, iron sound 
waves in plasma, and vibrations in a nonlinear string. 

Conclusions 

In this article, we have successfully formulated solitary 
waves solutions from the traveling wave solutions to 
the fourth order Boussinesq equation through the exp 
(-0(^))-expansion method. The procedure is simple, 
direct and constructive with the help of a computer 
algebra system. The method is quite efficient and prac- 
tically well suited to be used in finding solitary wave 
solutions of NLEEs and the attained solutions demon- 
strated the competence of the exp(-0(^))-expansion 
method. We also observed that the method is straight- 
forward and can be applied to many other nonlinear 
evolution equations. 
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